Thermal metastabilities in the solar core 
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Abstract 

Linear stability analysis indicates that solar core is thermally stable for 
infinitesimal internal perturbations. For the first time, thermal metastabil- 
ities are found in the solar core when outer perturbations with significant 
amplitude are present. The obtained results show that hot bubbles gener- 
ated by outer perturbations may travel a significant distance in the body 
of the Sun. These deep-origin hot bubbles have mass, energy, and chemical 
composition that may be related to solar flares. The results obtained may 
have remarkable relations to activity cycles in planets like Jupiter and also 
in extrasolar planetary systems. 
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"The Sun's interior is believed to be in a quiescent state and therefore the 
relevant physics is simple. However, after decades of work, we are still not cer- 
tain whether the physics of solar energy generation is as simple as expected" 
- wrote Bahcall |J. The linear analysis [Q, |3|, |J] indicated that solar core is 
stable to internal radial and non-radial perturbations. At the same time, the 
Sun is thought to be unstable for a wide range of magnetic and rotational 
instabilities and changing in a wide variety of timescales ||. Above a slightly 
larger than actual central temperature of the Sun, at 1.74 x 10 7 K CNO would 
begin to dominate over pp cycles. One may expect that outer perturbations 
generating local heating may set up convection-like movements in the solar 
core. Grandpierre |J suggested that planetary influences may interact with 
local magnetic fields in the solar core, and the arising local heating may gen- 
erate temperatures above 10 8 K. The possible role of solar inertial motions 
(SIM) around the mass centre of the solar system in the solar activity was 
recently emphasised JT], |8|, [|. In a frame fixed to the Sun the waves gener- 
ated by SIM have a relative velocity 1.5 x 10 5 cm/s, which is the rotational 
velocity of the solar surface. Their dissipated energies are derived from the 
rotational energy of the Sun. Solar rotation during the last 4.6 x 10 9 years 
has been spun down from the earlier 10 times the present rotation rate with 
Erot,o ~ 10 45 ergs, representing L ro t,o ~2x 10~ 6 Ls un - Although the present 
spun-down Sun has L rot ~ 2 x lCr 8 Ls un , it is still possible that at "elementary 
dissipation events" localised in space and time local rotation energy dissipa- 
tions may become significant at the most favourite sites where magnetic and 
planetary perturbation inhomogenities are the most pronounced. Consider- 
ing bubbles with VJ, > 10 18 cm 3 (i.e. the minimal bubble volume to reach 
solar surface, V c is the volume of the solar core ~ 10 31 cm 3 ), V c /Vb < 10 13 , 
the dissipation luminosity has an upper limit Ldi ss < L rot x 10 13 ~ 10 5 Ls un - 
Moreover, the episodic character of dissipation events offer a concentration 
in time in a rate C t ~ 1 to 10 6 leading to L^s < 10 11 Ls un . For a significant 
local heating to occur, we can allow to go down many orders of magnitude 
from the above derived upper limit. For example, with L^ss ~ 10 _1 £,swi, it 
is possible to reach a bubble initial temperature To > 1.7 x 10 7 K. 

The motion of a heated bubble is determined by the equality of the buoy- 
ant (F b = VgAp) and frictional (Ff = K/2v 2 Sp s ) forces, where S is the 
cross section of the bubble, V is its volume, K is the coefficient of turbulent 
viscosity, p s and Ap are the density of the surroundings and the density dif- 
ference between the bubble and its surroundings, and g is the gravitational 
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acceleration. Equating these forces, v 2 K/2S/V = gAp/p s . Now assuming 
pressure equilibrium between the bubble (referred with no index) and its sur- 
roundings (referred with index s), pT = p s T s , Ap/p s = (1 — T s /T). Taking 
K = l,we obtain for the bubble's velocity v = (8/3%(l - T/T s )f/ 2 . With 
typical values in the solar core g = 2 x 10 5 cms~ 2 , R = 10 5 — 10 6 cm, v ~ 2 
to 7 x 10 5 cms~ l jT(J. The estimation of thermal adjustment time is given 
by the following formula of [O] : 



Tadj = C p K P z R 2 /(AacT A ). (1) 

With typical values (« = 2cm 2 g~ 1 ,p = 90gcm~ 3 ,T = 10 8 K,R = 10 6 cm), 
Tadj = 3 x 10 3 s, while for T = 10 7 K, r adj = 3 x 10 6 s. We may define a 
surfacing time-scale for the bubbles 

Tsurf = h/v, (2) 

where It is the temperature scale height in the solar core It ~ 1.5 x 10 10 cm. 
With v = 1.5 x 10 5 cms _1 , r snr/ w 10 5 s. 

Estimating the nuclear heating time-scale of a hot bubble from e = 
e pT u = C p dT/dt, with a substitution x = T u , {l/x)dx/dt = eu/C p T, we 
obtained a relation |J to the local time-scale of the nuclear heating of the 
hot bubble as 

Tnud = C p T/eu. (3) 

When T changes from 10 7 K to 10 9 K, r nuc i changes from 10 16 s to 10 _1 s. We 
calculated the adiabatic expansion of the hot bubbles following Gorbatsky 
1(J, and the related time-scale is 



r ad = -5/(l/pdp/dr)l/(8/3(l - T s /T)Rgf' 2 (4) 

T a d is usually in the range of 10 4 to 10 5 s. From these estimations we can see 
that there are four relevant time-scales and these may be comparable for ther- 
mal perturbations in the solar core. Therefore, they indicate the possibility 
of thermal metastability and the need for detailed numerical calculations. 
We regarded that the material heated by the heat wave of radiative dif- 



fusion expanding from the bubble is coupled to it |10fl . We worked with a 
fourth order Runge-Kutta method to solve the differential equation system. 
They are based on: 

dR/dt = (-l/5)(l/pdp/dr)vR = R/r adj (5) 
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dQ/dt = (dQ/dt)ad + (dQ/dt) eps + em 



(6) 



dQ/dt « 2/hQ/pdp/dr + 3/2R g / pT s 4:TTpsR 3 /r adj + e4n/3R 3 p (7) 



8/3Rg(l -T s /T)f/ 2 (8) 



We used the opacity-figure of [12] and approximated it as k = io 10 - 45 /T 1 ' 45 
above 3 x 10 5 K. Our calculations differ from the previous ones that we con- 
sidered perturbations for which p' ^ 0, and T' ^ 0. This is a difference to 
the case considered elsewhere || |3|, |J with p' = 0, and (T/p)' = 0. In our 
calculations, we allowed the transient development of pressure inequilibrium. 
The duration of this pressure inequilibrium is a small fraction of a second, 
and therefore practically we have pressure equilibrium, p' = 0, (Tp/p)' = 0, 
but p' ^ and X" ^ 0. We solved the differential equation system with a 
numerical code developed by one of us (G. A.). We neglected the radiation 
pressure in all the terms except the diffusive one. Our method works well 
below and around 10 s K, the estimated errors in each quantities are smaller 
than 15%. At higher temperatures, the approach applied becomes more un- 
certain and therefore we could not follow the development of the hot bubbles 
at thermal runaways where our method indicated temperatures above 10 10 
K. 

We may observe in TABLE I. that above 10 8 K the heating time-scale 
becomes lower than that of the combined (expansion plus diffusion) cooling 
one. Our detailed calculations has shown that: 

1. ) Hot bubbles formed in the solar core may travel significant distances. 
Already with T = 1.74 x 10 7 K the bubbles take a path of w 1.2 x 10 10 
cm, i.e. even a moderate local heating may trigger hot bubbles travelling 
far from their site of birth. Depending on the size, energy and frequency 
distribution of initial perturbations, the solar core may become dynamic 
(meta-convection) . 

2. ) Hot bubbles may reach the solar surface even with Q = 10 35 ergs. 
The internal energy excess of the bubbles when arriving to subphotospheric 
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regions may reach 10 32 - 10 36 ergs, where the bubble's velocity is typically 
around 5-15 kms^ 1 . 

3. ) For initial temperatures higher than 10 s K numerical calculations 
indicate the presence of thermonuclear instability. 

4. ) Hot bubbles have a tendency to preserve their chemical composition 
corresponding to deeper layers, and consume a part of their hydrogen content, 
and therefore they will have "anomalous" chemical composition. 

Our results demonstrate the possibility of a coupling between some high 
velocity hot bubbles and some activity centres like rigidly rotating |L3| sunspot 
nests [|14], |15]1 , hot spots [nj, active longitudes |l6|j , and perhaps even with 
some individual atmospheric eruptions like flares, coronal mass ejections of 
protuberances. The possibility of such relations raises the question of relation 
between mass, energy and chemical composition of the bubbles and such at- 
mospheric events. The initial mass of a hot bubble is around 10 18 - 10 20 gram. 
Therefore, if its mass is conserved during its path to the envelope, in order 
to find a connection to flares with mass 10 15 - 10 16 g, a mechanism should 
exist between the suphotospheric and coronal regions to concentrate energy 
into a tiny fraction of mass. One candidate for such a mechanism is offered 
by the interesting fact that the bubbles approach the speed of sound when 
entering the subphotospheric regions. The developing sonic boom generates 
shock waves, transforming the energy of high speed bubbles into energy of 
shock waves and energetic particle beams [[L7j 



The bubble's cores when approaching the surface have a size of ~ 10 7 cm 
at r w 0.93-Rs„ n . In the subphotospheric regions they grow to ~ 10 8 cm, 
in good agreement with the photospheric sizes of flare kernels ~ 2000 km. 
Besides the energetic, mass and size relations, some connections between the 
solar neutrino problems, as well as the chemical anomalies of many flares 
and that of the hot bubbles may also be related. One may speculate that 
SIM inertial tides play an important role in the generation of solar activ- 
ity. The periodic acceleration and deceleration of SIM may drive an inertial 
oscillation in the solar core. Tidal stresses are maximal at the solar cen- 



ter [IB]. Magnetic fields may act as a restoring force, therefore magnetic 
torsional oscillation may arise which may explain the solar cycle ||. The 
maximum value of the Sun's orbital angular momentum reaches 25% of the 
Sun's rotational momentum 1.7 x 10 47 gem 2 s -1 . Spin-orbit coupling, a trans- 
fer of angular momentum from the Sun's orbit to the Sun's spin rate and 
vice versa could make a difference of more than 5% in the equatorial velocity 
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which is actually observed p 



Earth as well, they also show deep origin from the core p2 
chemical anomalies H 



Interestingly, hot spots occur at the 

rigid rotation 



|] and surface activity. A significant positive cor- 
relation exists between the chemically 'anomalous' activity of Jupiter's hot 
spots and the solar cycle ||25|| . The phenomenon of active longitudes is known 
also in close binaries where the position of the two diametrically opposing 
active longitudes is directly related to the position of the companion star 
(Olah, K. et al., in preparation). There is evidence for a purely tidal effect 
as well-enhanced activity at the subbinary point of close binaries f26fl . Tidal 
response calculations are also of interest in connection with the newly dis- 
covered planets J27j, which may play a role in the generation of the activity 
cycle of their host stars. 
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Table 1: The combined (expansion plus diffusion) cooling time-scale r coo ; in 
seconds for hot bubbles with different size and initial temperatures. For com- 
parison, the heating (nuclear) time-scale r nuc i is indicated at each tempera- 
ture below the relevant cooling times in the second column. T s = 1.318 x 10 7 
K is the temperature of the surroundings at r = 0.1Rs un from where the 
bubble starts rising in our calculations. Letter s indicates that the bubble 
reaches the solar surface, and letter r indicates the case when even thermal 
runaway occurs. 
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Figure 1: The run of the relative temperature of the hot bubbles for T = 
1.74 x 10 7 K, 5 x 10 7 K, 9 x 10 7 K and 2 x 10 8 K, when R = 3 x 10 8 cm. T and 
T(SSM) are the temperatures of the hot bubbles, and their standard solar 
model surroundings, respectively, r is the distance from the solar center. 
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Figure 2: The run of the hot bubbles' rising speed. The small waves around 
the end of our curves do not represent a phbysical effect as they are due to 
the linear interpolation applied in the solar model. 
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Figure 3: The run of the relative hydrogen abundances X/X(SSM) of the 
hot bubbles. 
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